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I.  INTRODUCTION 


I 

I 


Based  on  a fairly  extensive  set  of  empirical  data,  two  motion  sickness 
incidence  (MSI)  predictive  models  [A, 5]  have  been  developed  for  subjects 
exposed  to  sinusoidal  motion.  (In  essence,  MSI  is  the  probability  of  emesis.) 
One  model  [5],  which  will  be  referred  to  as  MSI  Model  I,  was  developed  for 
a fixed  exposure  time  of  two  hours.  The  other  model  [A],  which  will  be 
referred  to  as  MSI  Model  II,  includes  exposure  time  as  an  Independent  vari- 
able. Both  of  these  models  appear  to  yield  relatively  accurate  MSI  predictions 
in  those  cases  where  motion  is  defined  by  a single  sinusoid.  However,  neither 
model  can  be  used  directly  in  predicting  MSI  for  ship  motion,  because  that 
motion  tends  to  be  broadband.  Needless  to  say,  the  problem  of  obtaining 
accurate  MSI  predictions  tpr  actual  or  simulated  ship  motion  must  be  con- 
sidered. 

This  technical  report  summarizes  three  approaches  for  obtaining  such 
predictions,  and  examines  them  in  light  of  the  observed  MSI  data  from  an 
experiment  [3]  in  which  subjects  were  exposed  for  two  hours  to  motion  pro- 
duced by  the  sum  of  two  sinusoids.  These  approaches  will  be  referred  to  as 
(1)  the  independent  effects  approach,  (2)  the  Donnelly  weighting  approach, 
and  (3)  the  least  squares  weighting  approach.  Before  these  approaches  are 
defined  and  examined,  the  two  MSI  models  will  be  discussed,  and  the  experi- 
mental results  will  be  presented. 
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II.  THE  MSI  MODELS 


Because  the  two  MSI  models  [A, 5]  as  originally  stated  define  an  MSI 
value  as  a percentage,  a multiplicative  factor  of  100  is  required.  In  this 
report,  however,  the  alternative  (and  more  convenient)  expression  of  an  MSI 
value  as  a decimal  fraction  is  used.  Thus,  an  MSI  value  of  .426  in  this 
report  is  equivalent  to  an  MSI  value  of  42.6  (which  must  be  regarded  as  a 
percentage)  derived  from  an  MSI  model  in  its  original  form. 

Both  of  the  MSI  models  Involve  the  cumulative  distribution  function  $ 
of  a standard  normal  random  variable,  where 


-1/2  z 2 

'I'(z)  = (2ir)  f exp[-t  /2]dt 


MSI  Model  I [5]  predicts  MSI  for  two  hour  exposure  to  sinusoidal  motion  with 
frequency  f and  acceleration  a,  which  is  defined  as  "...  the  time  integral 
of  the  absolute  value  of  acceleration  imparted  in  each  half-wave  cycle."  The 
MSI  prediction  is  given  as 


MSI  = $[(w  - y )/a  ] 
w w 


where 


w = log  a 


y = y (f)  = .654  + 3.697  log  f + 2.320(log  f) 
w w 


and  0 = . 40  . 

w 


In  reference  [5],  the  value  of  a is  specified  as  .40  log  a . 
evidently  a misprint,  with  the  correct  value  being  » .40  . 


This  is 


J 


Because  MSI  Model  II  is  expressed  in  terms  of  rms  acceleration,  which  is  in 
more  common  usage  than  a,  MSI  Model  I will  be  transformed  into  terms  of  rms 
acceleration  (denoted  by  a)  for  comparison  purposes.  Thus,  in  place  of  the 
variable  w = log  i,  a new  variable  x = log  a is  required. 

As  indicated  in  reference  [5],  the  equality  a = .901  a holds  for 
sinusoidal  motion.  Therefore,  it  follows  that 

w = log  a 

= log(.901  a) 

= log  a + log (.901) 

= X - .045 

Thus,  by  transforming  from  variable  w to  variable  x In  equation  1,  MSI 
Model  1 may  be  expressed  in  terms  of  rms  acceleration  (a)  as 

MSI  = $[(x  - M^)/a]  (2) 

where 


X = log  a 

U,  - U,(f)  = U (f)  + .045 
11  w 

= .699  + 3.697  log  f + 2.320(log  f)^ 

and  a * 0 = .40  . 

w 

MSI  Model  II  [4]  predicts  MSI  for  t minutes  of  exposure  to  that  sinusoidal 
motion.  Its  predictions  require  the  definition  of  the  following  quantities: 
z = (x  - \i2)lo 

Zj.  " (y  - 

^t  ’ ^^t  " 

-3- 


where 


X = log  a 

2 

^2  = ^ ■‘’  2.73(log  f) 

a = .47 
y “ log  t 
= 1.46 
= .76 

and  p = - . 7 5 . 

The  resulting  MSI  prediction  is: 

MSI  - <I>(z)$(zp  . (3) 


H 


I III.  THREE  APPROACHES 

I 

Assuming  that  the  two  models  described  in  the  previous  section  provide 
accurate  MSI  predictions  for  motion  defined  by  a single  sinusoid,  a major 
question  is  whether  these  models  may  be  used  in  some  way  to  yield  MSI  values 
for  broadband  motion.  The  following  sections  outline  three  possible  approaches 
for  obtaining  such  predictions  from  the  MSI  models. 


A.  THE  INDEPENDENT  EFFECTS  APPROACH 

The  Independent  effects  approach  assumes  that  when  motion  results  from 
a combination  of  sinusoids,  there  is  no  interaction  (i.e.,  synergism)  between 
them  relative  to  MSI.  Let  S^  denote  a sinusoid  with  frequency  f^  and  rms 
acceleration  a^  . (This  notation  will  be  used  throughout  the  remainder  of 
this  technical  report.)  Under  the  assumption  of  independent  effects,  it  follows 
in  the  simplest  case  of  two  sinusoids  that 

MSI(S^,S^)  = MSI(S^)  + [1  - MSI(.S^)]  MSKS^) 

= MSI(S^)  + MSKS2)  - MSI(S^)MSI(S2)  (4) 

where  MSI(S^)  denotes  the  MSI  for  the  single  sine  wave  S^  and  MSI(S^,Sj) 
denotes  the  MSI  for  a combination  of  the  two  sine  waves  S^  and  Sj  . 

B.  THE  DONNELLY  WEIGHTING  APPROACH 

Donnelly  [1]  suggested  that  a qualitative  assessment  of  MSI  for  broad- 


frequency  f . He  proposed  defining  f to  be  the  "most  critical  frequency”, 
o o 

that  Is,  the  frequency  which  results  in  the  largest  MSI  prediction.  Based 
on  this  approach,  a weighting  function  w^  may  be  defined  such  that 

MS1(S  ) = MSI(S  ) 
i o 

where  a = a.w  . 
o 1 i 

To  determine  a motion  sickness  rating  when  a number  of  frequencies  are 
present,  Donnelly  used  the  weighted  rms  acceleration 

2 2 1/2 

a = 

o i i 

and  the  frequency  f in  a single  sinusoid  MSI  model.  It  should  be  noted 
o 

that  he  stressed  that  this  motion  rating  should  be  considered  qualitative 
rather  than  quantitative.  However,  there  is  little  doubt  that  there  will 
often  be  the  temptation  to  use  his  weighting  procedure  as  a quantitative 
measure  of  MSI. 

Although  Donnelly  emphasized  this  weighting  approach  for  data  expressed 
in  1/3  octave  format,  there  is  nothing  in  his  derivation  that  restricts  it 
to  this  situation.  For  motion  which  is  defined  by  the  sum  of  two  sinusoids, 
this  approach,  used  quantitatively,  would  define 


MSI(Sj^,S2)  = MSI(Sq) 


where  f is  the  most  critical  frequency 
o 

,2  2 2 2,  1/2 

a^  = (w^ai  + W2a2> 


The  weights  used  by  Donnelly  were 


log  w^  = -2.73(.77  + log  f^)‘ 


4- 


■ 


with  f defined  to  correspond  to  the  weight  w = 1.00  (i.e.,  f = .170  Hz), 
o o o 

These  weights,  however,  were  evidently  based  on  a partially  developed  version 


of  MSI  Model  II. 

For  the  two  models  considered  in  this  technical  report,  the  weights  cor- 
responding to  each  must  be  defined.  For  MSI  Model  I (equation  2), 

hi  = hi(fi) 

2 

= .699  + 3.697  log  f^  + 2.320(log  f^) 

= -.774  + 2.320C.797  + log  f^)^ 
from  which  it  follows  that 


log  w = -2.320(.797  + log  f^) 


(6) 


and 


f = .160  Hz  . 
o 


For  MSI  Model  II  (equation  3), 


^2  = 


.87  + 4.36  log  f^  + 2.73(log  f^) 


= -.871  + 2.730(.799  + log  f ) 

i 


from  which  it  follows  that 


log  w = -2.730(.799  + log  f^) 


(7) 


and 


f = .159  Hz  . 
o 


-7- 


E j nmnrrr  - 


C.  THE  EMPIRICAL  LEAsi'  SQUARES  WEIGHTING  APPROACH 


The  third  approach,  like  the  one  described  previously,  also  weights  all 

frequencies  back  to  a single  frequency  f . However,  whereas  the  previous 

o 

approach  specifies  the  weights  w^  a priori  based  on  an  existing  sinusoid 

model,  the  empirical  least  squares  weighting  approach  determines  them 

a posteriori  from  experimental  MSI  data  resulting  from  motion  corresponding 

to  complex  waveforms.  Thus,  this  approach,  unlike  the  other  two,  requires 

access  to  MSI  data  from  motion  that  is  other  than  sinusoidal.  The  following 

paragraph  summarizes  the  empirical  least  squares  weighting  approach. 

Given  motion  described  by  k frequencies  fj^,...,fj^  and  corresponding 

rms  accelerations  a^,...,aj^  , determine  (from  the  sinusoid  MSI  model)  the 

value  of  a that  corresponds  to  the  observed  MSI  for  frequency  f . Then  set 
° o 

2 2 2 2 2 

w,a,  + . . . + w a = a (8) 

11  k k o 

By  considering  n ^ k of  these  cases  (for  the  same  k frequencies),  values  of 

the  weights  may  be  determined  by  solving  the  resulting  set  of  equations  by 

least  squares.  Thus,  if  an  MSI  estimate  is  required  for  motion  described 

by  the  same  k frequencies,  but  not  necessarily  the  same  accelerations,  the 

values  of  w may  be  used  to  predict  (via  equation  8)  the  equivalent  a 
i o 

A predicted  MSI  value  may  then  be  obtained  by  substituting  a^  and  f^  in  the 
appropriate  MSI  model. 


t 


IV.  EXPERIMENTAL  RESULTS  AND  PREDICTIONS 


I 

i 

Table  1 sununarizes  the  results  of  an  experiment  [3]  involving  motion 
produced  by  the  sum  of  two  sine  waves.  In  addition  to  providing  observed 
MSI  data,  this  table  lists  the  frequencies  and  rms  accelerations  which  were 
reported  in  this  experiment.  (It  will  be  noted  that  results  are  given  for 
conditions  II,  III,  IV  and  V only,  since  condition  1 corresponded  to  a single 
sinusoid  used  as  a control.)  These  reported  frequencies  and  accelerations 
are  the  ones  used  in  this  technical  report  in  the  examination  of  the  three 
approaches  for  obtaining  MSI  predictions  for  motion  produced  by  complex 
waveforms. 

Table  2 summarizes  the  MSI  predictions  obtained  by  applying  each  of  the 
three  approaches  to  the  frequency  and  acceleration  data  of  the  four  conditions 
in  Table  1.  As  an  example  of  how  the  entries  in  this  table  were  derived, 
the  calculation  of  predicted  MSI  values  for  condition  II  based  on  MSI  Model  I 
is  outlined  in  the  following  paragraphs.  MSI  predictions  for  the  other 
conditions  and/or  for  MSI  Model  II  (using  t = 120  minutes)  may  be  obtained 
in  a similar  manner. 

Condition  II  is  defined  by  motion  which  is  the  sum  of  two  sine  waves 
Sj^  and  S2  . The  corresponding  frequencies  and  rms  accelerations  are 


Condition 


Table  1:  Results  of  the  Two  Sine  Wave  Experiment 
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Table  2:  Predicted  MSI  Results 


X = log (.14) 

= -.854 

M » p^(.16) 

2 

= .699  + 3.697  log(.16)  + 2.320[log(. 16) ] 

= -.774  . 

so  that 

MSI(S^)  - <!>[(-. 854  + .774)/. 40] 

= <!>(-.  200) 

- .421 

Similarly, 

MSKS2)  - <I>[(-.824  + .562)/. 40] 

= $(-.655) 

=»  .256 

Therefore,  from  the  independent  effects  approach  used  in  conjunction  with 
Model  I,  the  predicted  MSI  (from  equation  4)  is 

MSI(Sj^,S2)  = .421  + .256  - (.421)(.256) 

= .569 

Based  on  the  a priori  weighting  (equation  6)  for  MSI  Model  I where 

f » .160  , it  follows  that 
o 

log  w^  - -2.320[.797  + log(.16)] 

- -2.320(.797  - .796)2 


2 


log  = -2.320[.797  + log(.32)J 


= -2.320(.797  - .495)' 


= -.2116 


From  this  It  follows  that 


= 1.000 


w » .614 
2 


22  22  1/2 
a - [(1.000)  (.14)  + (.614)  (.15)  ] 

o 


The  predicted  MSI  (equation  5)  from  the  Donnelly  weighting  approach  Is  calcu- 
lated by  using  MSI  Model  I (equation  2)  where 


X “ log  a 


log (• 168) 


-.775 


U,  - M,(f  ) 
1 1 o 


U^(.160) 


-.774 


so  that 


r 


MSI(S  ,S  ) - MSI(S  ) 

12  o 

= $[(-.775  + .774)7.40] 

- $(-.003) 

- .498  . 

In  order  to  apply  the  least  squares  weighting  approach,  a set  of  observed 
MSI  results  corresponding  a number  of  cases  of  motion  described  by  the  same 
k frequencies  Is  required.  Since  each  of  conditions  II,  III,  and  V correspond 
(approximately)  to  frequencies  f = .16  and  ^2  “ * this  approach  may  be 

applied  to  the  data  from  these  three  conditions. 

I 

The  observed  MSI  for  these  three  conditions  were  .500,  .645,  and  .781  , 

respectively,  for  which  the  corresponding  values  of  a must  be  determined. 

o 

From  MSI  Model  I (equation  2),  the  value  of  a corresponding  to  an  observed 

o 

MSI  value  MSI^  Is  given  by 

$[ (x  - y)/a]  * MSI 


where 


X = log  a 


y = Ui(f„) 


y^(.l60) 


The  resulting  solutions  are: 


a - .168  for  MSI  = .500 
o o 


a = .237  for  MSI  = .645 
o o 


a - .344  for  MSI  = .781 
o o 


Therefore,  based  on  the  least  squares  weighting  approach  (equation  8), 


the  following  set  of  equations  results: 


wj(.14)^  +W2(.15)^  » (.168)^ 

2 2 2 2 2 
w^(.14)^  + W2<.14)  =(.237) 

w^(.08)^  + W2(.26)^  = (.344)2 


The  least  squares  solution  yields  the  weights 


w 


1 


(.325) 


1/2 


and 

w = (1.708)^^^ 

2 

By  using  these  estimated  weights  for  fj^  and  f2  , the  predicted  value  of 

= . 160  is  given  by 


be  used  in  MSI  Model  I (equation  2)  to 


1 

For  comparison,  the  weights  obtained  using  MSI  Model  II  are  w^^  • (.153 
1/2 

w,  - (1.548)  . 


a to  be  used  in  conjunction  with  f 
o o 


,22  2 21/2 

(wiSi  + W2S2) 


These  values  of  a and  f would  then 


obtain  predictions  of  MSI.  For  example,  based  on  these  weights  the  predicted 
MSI  for  condition  II  would  correspond  to 


j 


a^,  - t.325(.14)2  + i. 708(. 

- .212 

Thus, 

X » log  a 

o 

* log (.212) 

= -.674 

= M^(.160) 

- -.774 

so  that 

MSI(S^,S2)  = <t[(-.674  + .774)/. 40] 
- 4>(.250) 

= .599 
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V.  A STATISTICAL  EXAMINATION  OF  THE  MSI  PREDICTIONS 


In  this  section  the  observed  data  will  be  used  to  examine  the  validity 

of  the  three  approaches  which  have  been  discussed.  Under  the  hypothesis  that 

a given  approach  results  In  the  prediction  of  the  true  MSI  value  (denoted  by 

MSI  ) for  a duration  of  two  hours  under  specific  motion  conditions, 
o 

X - n(MSI  ) 

z = 

/ n(MSI^) (1  - MSI^) 

would  be  approximately  distributed  as  a standard  normal  random  variable,  where 
n » the  number  of  subjects  observed  under  this  motion  condlton 

and 

X “ the  number  of  subjects  who  were  overcome  by  emesis  within 

two  hours . : 

Based  on  the  observed  values  of  z,  a test  of  the  hypothesis  may  be  made 
using  the  data  summarized  In  Table  1 and  Table  2.  For  example,  the  results^ 
for  the  Independent  effects  assumption  using  MSI  Model  I are: 

16  - 32(.569) 

For  condition  II,  z = » -.788  , p * .431 

/ 32(.569)(.431) 

20  - 3K.547) 

For  condition  III,  z - » 1.098  , p ■ .272 

/ 31(.547)(.453) 

I 

- 

It  should  be  noted  that  since  a two-sided  test  is  appropriate  here,  the 
observed  p-values  reflect  this. 

-17- 
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For  condition  IV,  z = 


21  - 3K.485) 


= 2.144  , p = .032 


31(.485)(.515) 


25  - 32(.572) 

For  condition  V,  z = = 2.392  , p = .017 

/ 32(.572)(.428) 


Thus,  the  hypothesis  is  rejected  for  condition  IV  and  V at  a 5%  significance 
level.  It  should  be  noted,  however,  that  when  a number  of  tests  are  made, 
the  problem  of  multi-test  bias  arises.  Multi-test  bias  refers  to  the 
phenomenon  that  as  the  number  of  tests  made  increases,  the  chances  of  making 
at  least  one  type  I error  also  increases.  To  avoid  this  sort  of  contamination, 
it  is  desirable  to  use  a single  overall  test. 

Such  a test,  developed  by  Fisher  [2],  may  be  based  on  the  statistic 


X = -2  E ln(p  ) 
i=l  ^ 


where  p denotes  the  probability  of  the  observed  outcome  of  the  1 — individual 
test  in  a total  of  N tests.  Under  the  hypothesis,  the  statistic  X has  an 
approximate  Chi-square  distribution  with  2N  degreees  of  freedom. 

For  the  data  under  consideration, 

X = -2[ln(.431)  + ln(.272)  + ln(.032)  + ln(.017)] 

= 19.348  (p  = .014) 

Thus,  there  is  convincing  evidence  that  the  independent  effects  approach,  based 
on  the  single  sine  wave  MSI  Model  I,  cannot  be  used  to  predict  MSI  for  a com- 
bination of  sine  waves.  It  follows,  therefore,  that  either  (a)  the  single 
sine  wave  MSI  model  is  Incorrect,  (b)  the  independent  effects  approach  is 
Incorrect,  or  (c)  both  are  Incorrect. 


• ■ * 


Tests  involving  the  Donnelly  weighting  approach  and/or  MSI  Model  II  may 
be  carried  out  in  a similar  manner.  Table  3 summarizes  the  results  of  these 
tests.  As  can  be  seen  from  this  table,  the  experimental  results  shed  serious 
doubt  on  the  usefulness  of  both  the  Independent  effects  approach  and  the 
Donnelly  weighting  approach  as  methods  of  predicting  MSI.  Tests  of  the  least 
squares  weighting  approach  are  inappropriate  and  therefore  not  Included  in 
the  table.  For  that  approach,  extremely  good  agreement  between  observed  and 
predicted  MSI  values  is  to  be  expected  because  two  weights  were  estimated 
based  on  a total  of  only  three  observations.  Furthermore,  predicted  MSI 
values  would  be  compared  with  the  observed  MSI  values  from  which  the  weights 
had  been  estimated.  A valid  test  of  this  approach  would  require  a larger  set 
of  observed  data. 

It  is  interesting  to  note,  however,  that  the  weights  estimatevi  by  this 
approach  are  reversed  in  magnitude  from  the  weights  determined  by  the  approach 
suggested  by  Donnelly.  Whereas  the  weight  his  approach  assigns  to  a frequency 
is  a monotonic  decreasing  function  of  the  distance  of  that  frequency  from 
approximately  .160  Hz,  the  weight  assigned  to  .330  Hz  by  the  least  squares 
approach  is  greater^  than  the  weight  assigned  to  .160  Hz. 

Not  enough  data  exists  with  which  to  test  the  use  of  MSI  predictions  from 
the  Donnelly  weighting  approach  as  qualitative  (i.e.,  ranking),  rather  than 
quantitative,  indices.  However,  the  observed  and  predicted  MSI  values  listed 
in  Table  2 indicate  that  the  Donnelly  approach  ranks  (in  order  of  increasing 
MSI)  conditions  IV,  III,  II,  V based  on  Model  I,  and  conditions  IV,  III,  V,  II 
based  on  Model  II,  while  the  observed  results  indicate  the  ranking  II,  III,  IV, 
V,  Instead. 


1 

Two  times  greater  using  MSI  Model  I and  over  three  times  greater  using 
MSI  Model  II. 
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Condition 

Independent  Effects 
Approach* 

Donnelly  Weighting 
Approach* 

II 

z = -,788(-1.155) 

z = .023(-.317) 

p = .431(  .248) 

p = .982(  .751) 

III 

z = 1.098(.768) 

z = 1.784(1.450) 

p = .272(.442) 

p = .074(  .147) 

IV 

z = 2.144(2.121) 

z = 2.795(2.616) 

p = .032(  .034) 

p = .005(  .009) 

V 

z - 2.392(2.242) 

z = 3.058(2.980) 

p = .017(  .025) 

p = .002(  .003) 

Overall 

X = 19.348(18.570) 

X = 27.970(25.548) 

p = .014(  .018) 

p = .001(  .001) 

A 

The  first  value  Is  obtained  using  MSI  Model  I (equation  2).  The 
value  in  parentheses  is  obtained  using  MSI  Model  II  (equation  3). 


I 

Table  3:  Summary  of  Statistical  Tests 


-20- 


1 

I 

A 


VI.  SUMMARY  AND  DISCUSSION 


Although  only  a small  amount  of  empirical  MSI  data  exists  for  motion 
other  than  that  produced  by  a single  sinusoid,  that  data  provides  a strong 
indication  that  neither  the  Independent  effects  approach  nor  the  Donnelly 
weighting  approach  produces  accurate  MSI  prediction.  (Judgment  about  the 
least  squares  weighting  approach  must  be  reserved  until  enough  empirical 
data  exists  to  provide  an  adequate  test.) 

It  should  be  pointed  out  that  the  statistical  analysis  carried  out  in 
this  technical  report  was  approximate.  Although  the  analysis  did  use  the  pair 
of  observed  frequencies  and  their  associated  rms  accelerations  as  stated  in 
reference  [3],  use  of  the  actual  frequency  spectrum  might  have  produc  d slightly 
different  results.  Furthermore,  the  analysis  assumed  that  the  predictions 
obtained  from  MSI  Model  I or  II  were  true  MSI  values,  devoid  of  variability. 
Since  the  development  of  each  of  these  models  was  based  on  a relatively  large 
number  of  observations,  the  effects  of  such  variability  should  be  negligible. 

In  any  event,  with  results  as  extreme  as  those  observed,  it  is  highly  doubtful 
that  these  two  approximations  would  have  changed  the  conclusions  which  were 
reached . 
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This  report  summarizes  three  approaches  which  have  been  suggested  for 
predicting  motion  sickness  incidence  (MSI)  for  actual  or  simulated  broadband 
ship  motion.  Two  of  these  approaches  are  based  on  the  results  of  MSI  models 
developed  for  pure  .sinusoidal  motion.  The  third  approach,  which  uses  weights 
derived  by  least  squares,  cannot  be  based  directly  on  these  models.  Instead, 
it  would  be  developed  empirically  from  observed  MSI  in  experiments  involving 
broadband  motion.  — fcontlnued 
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20.  Abstract  (continued) 

Although  only  a small  amount  of  empirical  MSI  data  exists  for  motion  other 
than  that  produced  by  a single  sinusoid,  this  data  provides  strong  statistical 
evidence  that  the  two  approaches  based  on  the  single  sinusoid  models  fall  to 
produce  accurate  predictions.  Statistical  judgments  about  the  least  squares 
weighting  approach  must  be  reserved  until  enough  empirical  data  exists  to 
provide  an  adequate  test. 
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